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Abstract. We study dynamics of interfaces in solutions of the equation ^f^{u) — 0, 

for /e of the form /e(it) — (it^ — l)(2it — ek), for k G R, as weU as more general, but 
qualitatively similar, nonlinearities. We prove that for suitable initial data, solutions exhibit 
interfaces that sweep out timelike hypersurfaces of mean curvature proportional to k. In 
T-H particular, in 1 dimension these interfaces behave like a relativistic point particle subject 

to constant acceleration. 

(N 

1. INTRODUCTION 

(N 

In this paper we consider the dynamics of interfaces in semihnear hyperbohc equations. 
p I The simplest example that we study is the equation 

(1.1) e{uu-u^x) + -{u'^-l){2u-eK) = 0, 

where we assume for concreteness that k > 0. Here the nonlinearity /e(n) = (n^ — l){2u — eK) 
^ has the form = F^, where has local minima at u = ±1, with ^^(ibl) = ib|e«;. Thus the 

I— I state u = —1 has slightly lower energy than the state u = 1, and one might expect that there 

exist solutions in which the low-energy phase u = —1 grows at the expense of the higher- 
^ energy phase. This is what we prove. In fact we show that, for suitable initial data, solutions 

^\ exhibit an interface that behaves like a relativistic mass subject to constant acceleration 

proportional to the parameter k. Equivalently, the interface sweeps out a timelike curve of 
constant Minkowskian curvature, proportional to k, in the (t, a;)-plane. 

It turns out that our analysis extends with rather few changes to wave equations on 
suitable Lorentzian manifolds (A^, h). Thus we will also consider the equation 

cn 1 

^ (1.2) enhU + -/o(u) + K/i(u) = 

• • e 

where Dh is the Laplace-Beltrami (wave) operator on (A^, h), k is a smooth function, and 
^ ^foi'iJ') + i'^fi{u) generalizes the nonlinearity in (1.1) in a natural way. In this situation, 

^ analogous to ( |1.1[ ), we show that for well-prepared data, interfaces sweep out timelike hy- 

persurfaces of prescribed mean curvature k, with respect to the Lorentzian metric h. 

In the case when (A^, h) is just 1+n-dimensional Minkowski space and k = 0, corresponding 
to the situation when the two potential wells have equal depth, similar results were proved 
by the first author in following partial results of Thus, the present paper consists of 
a number of improvements of the basic argument developed in [9]: we extend the results to 
the case k 7^ 0, we show that they remain valid on Lorentzian manifolds more general than 
Minkowski space, and we drop some convenient but artificial restrictions imposed in [9j on 
the topological type of the hypersurfaces considered. A key point in our analysis is that if k is 
a nonzero constant, then in certain weighted energy estimates it is much more useful to use. 



O 



not the canonical conserved energy associated to the actual equation (1.2) under study, but 



rather the conserved energy associated to the k = Q equation. (See Remark 2.1 ) This simple 
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observation plays a crucial role in our arguments and makes the extension of techniques 
developed in [9] to the more general situation considered here surprisingly straightforward. 

Equations such as (1.2), with k ^ 0, have been studied in the cosmological literature as 
models for what is called the decay of a false vacuum. This arises from models in which the 
universe is described by a quantum field theory for which an equation like (1.2) (or a more 
complicated but in some ways similar equation) is a low-energy limit, and whose state is 
initially given by a constant function u = Vf, where vj is a "false vacuum": a local, but not 
global, minimum of some underlying potential function. In the example (1.1), if k > then 
Vf = 1, and the "true vacuum", or global minimizer of the potential function F^, is vt = —1. 
In this situation, a quantum tunnelling event could in principle lead to the nucleation of 
region in which u = vt- This scenario was investigated in a series of papers by Coleman 
and coworkers; see for example |i4J, which estimates via a formal semiclassical approximation 
the probability per unit time per unit volume of such a tunneling event. Our results have 
nothing to say about this, but describe the dynamics of a fully-formed interface between 
false and true vacuums in a universe governed by (1.2), showing that if the interface has an 
energetically optimal structure, then it behaves like a hypersurface of constant Lorentzian 
mean curvature proportional to the difference in energy between the true and false vacuums. 

Earlier work on dynamics of energy concentration sets in hyperbolic equations includes 
[HI \TT\ [6l [16], all of which consider situations in which energy concentrates around points 
rather than submanifolds, as in [31 [9] and the present paper. The dynamics of interfaces in 
equations such as (1.2) is studied from a formal point of view in [14j. 

A lengthy discussion of related elliptic and parabolic results, with a heavy bias toward 
the K = case, is contained in [9]. For k ^ 0, there is a rather strong analogy between the 
phenomena we study and propagating fronts in semilinear parabolic equations, a subject that 
has attracted a great deal of study, dating back to the 1930s [TO]. In particular, the problem 



that formally determines the profile and (relativistic) acceleration of interfaces (see (2.5) or 
( 3.17| )) is exactly the same one that determines the profile and velocity of propagating fronts 
in certain parabolic problems, see for example [HIS [2]. There is also an analogy between our 
work and results that establish an asymptotic connection between elliptic analogs of (1.2) 
and surfaces of prescribed Euclidean (or more generally Riemanian) mean curvature, see for 
example [13l E] . 

This paper is organized as follows: In order to highlight some main ideas with as few 
preliminaries as possible, we consider in Section [2] the case (1.1) of an equation in one space 
dimension associated with two potential wells of unequal depth. This discussion isincluded 
just to illustrate our arguments in a simple setting, and is not needed in later pars of the 
paper. 

We therefore defer until Section [3] both the statement of our main result, and the intro- 
duction of some notation that is used throughout the rest of the paper. In Section |4] we 
introduce a coordinate system in which many of our main estimates will take place, stating 
the properties that we will need and deferring most proofs to Section [7} The heart of our 
argument consists of weighted energy estimates in this adapted coordinate system. These 
are carried out in Section [5} In Section [6| these estimates are combined with rather standard 
energy estimates away from the interface in a iterative argument that completes the proof 
of our main theorem. 



2. THE SIMPLEST NONTRIVIAL EQUATION 



In this section we consider the 1-dimensional equation ( 1.1 ). All the results in this section 
are essentially subsumed in Theorem 3.2 and most of the main ideas in Theorem 3.2 appear 
here, in somewhat simpler form. 
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It is convenient to consider initial data such that at t = 0, 

(2.1) (u, ut) = (-1, 0) for all x near 0, [u, ut) = (1, 0) for all x > R 

for some R. (More conditions on the data will be imposed later.) Noting that the constant 
functions ±1 are both solutions of (1.1), standard facts about finite propagation speed for 
solutions of (1.1) imply that 

(2.2) w(t,a;) = -1 for (t,x) near (0,0), u{t, x) = I for x > R + \t\. 

2.1. change of variables. As suggested above, one might guess that for suitable initial 
data, solutions will exhibit an interface that sweeps out a timelike curve of constant (nonzero) 
Minkowskian curvature proportional to the parameter k that controls to the difference in 
depth of the two energy wells. Such curves have the form {{t,x) : x^ — = c,x > 0} 
modulo translations and reflections. We thus start by changing variables in such a way as to 
"straighten out" a 1-parameter family of such curves. Thus, we introduce new coordinates 
{e, r) defined (for 6* G M, r > 0) by 

(2.3) {t,x) = (rsinh6',rcosh6') = ip{9,r) £ { {t,x) : |t| < x}. 

These are just Minkowskian polar coordinates, with 6 being the angular and r the radial 
coordinate. Note that every coordinate line r = rg is a timelike curve of constant curvature 
^ with respect to the Minkowski metric, which in these coordinates takes the form ds"^ = 
—r'^dO'^ + dr^. We will treat ^ as a timelike coordinate, and r as spacelike. 
If u solves (1.1) and we write v = uotp^ then we find that v satisfies 

(2.4) e{^vee-Vrr--vA + ^{v'^-l){2v-eK) = Q, 6'GM,r>0. 

r J £ 

If we imagine that vqq ~ and that ^ ~ c constant, then this looks like the equation 

(2.5) £{-q" -cq')+^^{q^ -l){2q-£K) = 

for the profile q and wave speed c of traveling wave solutions of the parabolic counterpart of 
(1.1). This is known to have the 1-parameter family of solutions 

T — Tn 

(2.6) c = K, g = tanh(— -^), tq G M. 

(We have set things up so that the profile q is independent of the parameter k.) Note that if 
we choose = ^, then all the nontrivial behavior of qe{r) := tanh(^(r — ^)) is concentrated 
in an e- neighborhood of r = -, which is consistent with the heuristic - ~ k = c. 
Thus, we will study (2.4) with initial data such that 

(2.7) ?;(0,r) PS tanh(^(r - ro)), ^^(O, r) for r > 
where henceforth we set 

(2.8) ro := ^ 

K 

Indeed, we will show that for data of this form, solutions are approximately independent of 
6, and hence remain concentrated about the curve r = \. Recall also that we are assuming 
(2.1), which in the new variables implies that for ^ = 0, 

(2.9) {v, ve) = (-1, 0) for aU r close to 0, {v, ve) = (1, 0) for r > R, 



and this implies (2.2), which translates to 

(2.10) i;(6',r) = -1 for r near 0, v{e ,r) = I ioi r > Re^'^^ . 

for every G M. 
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2.2. differential energy inequality. We next define 

1 



I IE 



A short computation shows that if is a sufficiently smooth solution of (2.4), then 

d 



(2.11) 
where 



eeiv) = £{v0Vr)r + Term 1 + Term 2 



Term 1 = e—Vr(l — Kr 



Term 2 = Keve{vr (1 — v )). 



Term 1 should be small in if vq is small and Vr is concentrated near r = ^. Also, the 
profile qeir) = tanh(^(r — ^)) satisfies drQe = ^(1 — Q^), so that if f w in a sufficiently 
strong sense, then Term 2 should be small. 



Remark 2.1. Note that (2.4) has an exactly conserved energy: a computation shows that 
d 



dO 



{j^vl + ^ -^F,{v)^ = {rvrve)r, where F^{s) = (s^ - l){2s - ke) 



It turns out that it is much more useful to work with the approximately conserved energy 
ee{v) defined above. This observation, although very simple, is a key point in our analysis. 



2.3. lower energy bound. Note that if v satisfies (2.10), 

1 







lYdr > 



> 







|(l-t;^K| 

' 1 „ 
(v - -V )rdr 



(2.12) 



=: Co. 



2.4. weighted energy estimates in new variables. Next, given a solution v of (2.4), we 
will write 

\2l 



(2.13) 



cm 



[1 + (r - ro) ]ee{v) dr 







Co- 



For the initial data we consider, (2.10) holds, and then ( |2.12 ) implies that 
Ci(^) > (2(9) := r li-f + ir- ro? ( U + - 1 



(2.14) 



dr 



Using (2.11), we compute 

/>oo 

C[{e)= / [1 + (r - ro)^] [e(t;et'r)r + Term 1 + Term 2] dr. 
Jo 

Every term in the integrand contains a factor of vg, which due to (2.10) has compact support 
in (0, oo), so the integral clearly exists, and we can integrate by parts without problems. It 
also follows from (2.10) that 1 < Re^^^^ on the support of vg. Thus 



[l + ir - rof]E{vgVr)r 



2EVg{r -ro)vr 



< Re\'\ / 2e 



\ve\ 



rol l^r-l < 2i?el^lC2(e). 



Recalling that k = r^^, elementary estimates yield 



I Term 1| < 



EK 



2 



(^)^ + (r-ro)V 
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and 



iTerm 2| < 



rsK 



r 



1 



Repeatedly using (2.10) to bound r, and recalling that ^2 < Ci) we deduce that 

1 



ClW <Ci?e2|^ICi(0) + C7i?3e3|e| 



However, arguing as in ( 2.12[ ), 

1 



^ 2 , 2\2 
-V^ + -{l-V ) 



l |2.10| l 



^ 2 , 2\2 

-v^ + -{l-v ) 



Co 



(2.15) 
We conclude that 

(2.16) C'M < C{R + i?3)e3|''ICi(^). 

2.5. conclusions about v. At this point, we have proved most of the following proposition. 



Proposition 2.2. Let v solve (2.4) with initial data satisfying (2.9). Let vq = , and 

define Ci, C2 in (2.13), (2.14). Then there exists a constant C , depending on the parameter 
R in (2.9), such that for every G M, 

C2{e) < cm < e^(^"'"-i)Ci(0) for every 9, 



(2.17) 

As a result 
(2.18) 



dr 



j{9,r)-v{0,r)r ^ < 0(9, R) 



,Ci(0) 







In particular, there exists initial data for which the solution v satisfies 



Ci9,R):=C\9\ / e^(«)(^"^'-i)ds. 







(2.19) 



sign(r-ro)P^ < C(^,i?)e for all 9. 



Remark 2.3. For any 5 > 0, there exists C = C{5,ro) such that 

.r-ro . 2 dr 

. )-sign(r-ro) ^-^tk 

£ max(r, 0)^ 

ro 2 dr 



(2.20) 

So ( 2.19| ) implies that 
(2.21) 



|f (0, r) — tanh(- 



max(r, 5)^ 



< Ce. 



However, (2.20) says precisely that (2.19) is not a sharp enough estimate to determine the 
profile of v, so it would arguably be a little misleading to insist on tanh rather than sign in 
estimates such as (2.21), (2.19). 

On the other hand, standard spectral estimates imply that for every 9, there exists some 
ri;{9) such that the solution v in (2.19|) satisfies 



\v{9, r) — tanh( 



r — rir(9) ..2 dr 



< CCi{9) < Ci9)e' 



and then it follows from (2.19) that \rs{9)—ro\ < Ce. So although it is not captured in (2.19), 
our estimates do in fact show that v is close to a scaled, translated hyperbolic tangent. 
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Proof. To complete the proof, notice first that (2.17) follows from (2.16) and (2.14), via a 
form of Gronwall's inequality. Next, for every r. 



(v{9,r) -v{0,r)J < \e\ j -^vl{s,r)ds. 



We deduce (2.18) by integrating this inequality with respect to r then using Fubini's Theorem 



and (2.17). Finally, to prove (2.19), it now suffices to exhibit initial data {v, vq)\q=q satisfying 



(2.9), and such that 

Ci(0) < Ce\ 

0, and let 



(2.22) 
To do this, let ^0161=0 



/ |f (0, r) — sign(r — ro)|^^ < Ce^. 
Jo ^ 



w(0, r) : = 



where 

(2.23) qe,r,{s) = 

and for r > we define Xr £ C, 

(2.24) Xr 



Xro('S)tanh(-) + (1 



Xr-o('S))sign(s) 
) to be a function such that 



s) = 1 if |s| < I and Xr{s) = if \s\ > y. 



I I 



< C/r 



It is straightforward to verify that this initial data satisfies (2.22). 



□ 



2.6. conclusions about u. Proposition |2.2| yields uniform estimates for v on sets of the 
form {{r,9) : r > 0, |^| < ©}, corresponding to uniform estimates of the original solution u 
in a sector {{t,x) : x > 0, |t| < xtanh©}. (Recall that u = v o ip {oi ip defined in (2.3).) We 
next show infer estimates of u in a spacetime slab (— T, T) x M. 



Proposition 2.4. Fix e £ (0, 1] and let u solve (1.1) with initial data n(0,x) = (le,ro{x — vq) , 



Uj(0,2;) = 0, where q^^ra is defined in (2.23) above and r^ = k 



For every T > 0, there is a constant C{T), independent of e, such that 
t-T 

(2.25) / \u{t,x) - sign{x - -f{t))\^ dxdt < C{T)e. 



-T 



where ^(t) 



(r2 + t 



2\l/2 



Remark 2.5. Although we have stated here only an analog of (2.19), our arguments also 



establish an analog of (^2.17), i.e., energy estimates showing that, for a large class of initial 



data, energy concentrates near the curve (i,7(t)). We already know this, modulo a change 
of variables, in the sector {(i, x) : x > 0, \t\ < xtanh B} so the new point is energy estimates 



outside this sector, which are established in (2.28| below. 



Proof. Fix T > and let G be such that T = ^sinhG, i.e., 9 = sinh"^ ^. We will start 
by considering a solution u for initial data satisfying 

(2.26) (ii,nt)(0,x) = (-1,0) for all X < (n, uj)(0, x) = (1, 0) for all x > i?, 

for some < 6 < rQ < R. We will only specialize later to the initial data in the statement 
of the Proposition. Until we do so, all constants in our argument may depend on rg = ^ 
(hence on T), and the parameter R above. 

u o ijj^^ 



e 



The choice (2.26) of initial data implies that Proposition 

xi{t) : = 



2.2 



applies to v 



We will write 



xo(t) := \ t^ + 



9 



0' 



t2 + 



'0- 



The idea is simply that results from Proposition 2.2 imply that u ~ — 1 (with respect to 
the norm) in the set {(t, x) : xo{t) < x < xi(t)}, see the shaded region in Figure [l] We 
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combine this with the fact that {u, ut)\t=o ~ 0) to argue that u ~ —1 in in the entire 
set {{t,x) : |t| < T, X < xi{t)}. 




Figure 1. 



Step 1. Indeed, let x be a smooth function such that < x ^ 1 Suiid 
x{x, t) = 1 if a; < xo{t), xi^j t) = ii x > xi{t). 
We will write e°(n) for an energy density in the original coordinates defined by 



1 



Then rather standard energy arguments, which we recall in Step 3 below, imply that 

ft fXl{t) 



(2.27) 



Xe°(n) 



< e 



Kt 



+ Ci 



Jxo{t) 



3''(*-")e°(n) dx ds. 



It is straightforward to check that there exists some C such that e°(n) o < Ce^ {v) in 
[— 0, 0] X [|ro, .|ro]. On the same set, the Jacobian determinant det(-DV) is bounded and 



C'-y-' - 2' r 
Hence by a change of variables, 

ft fXl{t) 

e°Ju)dxds < C 



2£ 



Jxo{t) 



e /.2ro/3 



ee{v) dr de<c C2{0)de. 



Jro/3 







We combine this with (2.27) and using (2.17), and note that / xi^j x)e'^{u){0, x) dx < CCi(O), 
as a result of (2.1). These computations lead to the inequality 



(2.28) 



{{t,x):\t\<T,x<xo{t)} 



e°{u) dxdt < CCi(O). 



Step 2. We now use the above estimates to prove (2.25) for u solving (1.1) with initial 
ta si 

(2.29) 



data satisfying (2.26) together with 

Ci(0) < Ce^ ioxv = uoip~^, I |ti(0, x) — sign(x — ro)|^ < Ce 

JR 
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In particular, these conditions are satisfied by the specific data described in the statement 
of the proposition, exactly as in the proof of Proposition 2.2 We will write Z(t,y) := 
{t, y + xo{t)). Then for for every y, 



uiZ{t,y))-uiZiO,y)) = / utiZ{s,y))+x'ois)u,iZis,y)) ds < C ( t 



' et{u) {Z{s,y)) 



Also, {{t,x) : |t| < T, X < Xo{t)} = {Z{t,y) : \t\ < T,y < 0}, so by squaring the above 
inequality, integrating from y = —oo to y = 0, integrating in t, changing variables, and using 
( |2.28| ) and ( |2.29D , we find that 

f-T f-xo{t) n 
/ / \u{t,x) + l\^dx dt < -Ci(0)<Ce. 



(2.30) 



Next, recall that u o -0 = so by a change of variables (see (2.10)) and (2.19), 

rT j-R+\t\ 



-T Jxo(t) 



In — sign(x — 7(t))| <C 



rRe 



e Jro/3 



|f (^, r) — sign(r — ro)| dr d9 <Ce. 



Recalhng that u{t, x) = I ior x > R + \t\, we deduce that if ( |2.26D , ( |2.29D hold then 

rT 

|2 



(2.31) 



|u(t, x) — sign(x — 7(t))| <Ce. 



Step 3: proof of (2.27): A short calculation shows that 

d 



(2.32) 



dt 



e"{u) = e{utUx)x + nut{l -u). 



Then it follows from (2.32) that 



We integrate by parts and use some elementary inequalities to find 

d 



dt 



J xel{u) < J {\xt\ + |Xx|)e° + kJ xe°. 



This implies that 

d 
dt 



Kt 



Xeliu) < e-«* / {\xt\ + |Xx|)e°(n) < Cie" 



Kt 



for \t\ < T, where we may take Ci = — -. We arrive at (2.27) by integrating this 
expression from to t. 



□ 



3. STATEMENT OF MAIN THEOREM 

In this section we state our main theorem, which relates a semilinear wave equation to 
a hypersurface of prescribed mean curvature on a Lorentzian manifold. We first introduce 
these ingredients. 

We will always use greek letters such as a, /3, ... to denote indices that run between and 
n, and we implicitly sum repeated greek indices from to n. We will explicitly indicate sums 
that run over different ranges, we will generally not use greek letters for indices that belong 
to some proper subset of {0, . . . , n}. 
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3.1. the Lorentzian manifold. We consider a manifold that we assume to be home- 
omorphic to R-'^^", n > 2, and we fix global coordinates {x^, . . . ,x"). We will write h to 
denote a Lorentzian inner product on A^, and we let (ha/s) denote the components of the 
metric tensor with respect to the given coordinates. We also write 

= (/ia/3)~\ h = det{h^f,). 

We will assume that (ha/s) is smooth and that there exists some constant ci > such that 



(3.1) 



hoQ < -ci, 



< ^ for all a, /3. 



everywhere in N. Thus may be thought of as a time coordinate, and we will sometimes 
write t for x^. We will further assume that 



(3.2) 



hoi = 0, 



1, . . . ,n. 



Then it is clear that 0<c< —h<C everywhere in A^. For <t <T < oo, we will use the 
notation 

Et:={x£ N -.x^ = t}, NT:={xe N : \x°\ < T}. 



In view of (3.1), we can obtain a Riemanian metric from h by changing the sign of h 



00- 



Since this metric is uniformly equivalent to the Euclidean metric, and the associated volume 
element is uniformly comparable to the Euclidean volume element on M^+", we will for 
simplicity use the Euclidean structure on ]R^+" to define and Sobolev norms on N. 



3.2. the semilinear wave equation. Let k : N 

consider the equation 



(3.3) 
where 



0, 



be a fixed smooth function. We will 



u : N 



-h 



fiU 



It follows from (3.1) that (3.3) is a hyperbolic equation. We always assume that the nonlin- 
earities /o, /i in (3.3) have the form 

F\ h = 



(3.4) 

for F 
(3.5) 



/o 



a smooth function such that 




in [-1,1] 
elsewhere . 



F{x) > if |x| / 1, c(l - \x\f < F{x) < C(l - \x\f if \x\ < 2. 



Note that /o, /i are smooth as a consequence of ( |3.5| ) and the smoothness of F. 

We do not address questions about well-posedness of (3.3). When (N, h) is flat Minkowski 
space M}~^^ then global well-posedness in the energy space can be guaranteed by imposing 
suitable growth conditions on F. In particular, if n < 4 then (3.3) on Minkowski space is 
globally well-posed for fo{u) = 2{u^ — l)u and fi{u) = 1 — v?, associated to the potential 
F{u) = \{l-u^f. 

The form of the nonlinearity is further discussed in Section 3.6, where we show that 
the assumptions (3.4) are not actually restrictive if k is constant. Note also that one can 
easily generate nonlinearities satisfying the above conditions by starting from /i such that 
sign/i(s) = sign(l - s2), |/((±1)| / 0, then deflning F = \fl and h = F' = hf[. 
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3.3. the hypersurface of prescribed mean curvature. We assume that I is a bounded 
open subset of Nt* , for some T* > 0, such that T := dlDNT* is a smooth embedded timelike 
hypersurface satisfying the prescribed mean curvatur^ equation 

(3.6) — k{x) = mean curvature in {N, h) of F at x 
and that 

(3.7) r is orthogonal to the initial hypersurface Sq. 



This means that F has zero initial velocity with respect to the initial hypersurface. See (4.2 ) 



(4.3) below for a precise formulation. 

For smooth data and smooth k, such as we consider here, local existence of smooth 



embedded submanifolds F C satisfying (3.6), (3.7) follows from arguments in Milbredt 
Remark 3.1. In fact Milbredt [I2j studies a rather general Cauchy problem for the k = 



case of (3.6) on a larger class of Lorentzian manifolds than we consider here. Modifying 



his arguments to extend to the case of smooth k presents no difficulty. His basic existence 



results are ultimately proved by solving (3.6) in coordinate charts (with a suitable choice of 
gauge) and then piecing together these local solutions, using finite propagation speed for the 
equation. Solvability in coordinate charts depends on local solvability results for the Cauchy 
problem for a general quasilinear hyperbolic equation, of the form 

where g^^ < —X and g^^ > fi5^^ . Changing the equation to allow nonzero k simply adds some 
additional smooth lower-order terms on the right-hand side, and the existence results used in 
[T2j apply with no change to the equation once it is modified in this way. Other aspects of 
the argument, such as piecing together local solutions, are similarly uneffected. 

3.4. main theorem. The main result of this paper is the following. The statement uses 



terminology and notation introduced in Sections |3.1| - [3^ and for any set C/, we will write 

signjj{x) :-- 



1 if X G [/, 
— 1 if not. 



Theorem 3.2. Assume that Tq < T* . 

Then there exists a neighborhood N' ofTDNxQ and a unique smooth function d-p ■ N' ^ M 
such that 

(3.8) drix) = onT, h""^ d^c^dr d^pdr = 1, dr > in I n N'. 

Moreover, dr is bounded away from outside of every neighbourhood o/F. 

In addition, for every e G (0,1], there exists smooth initial data (uo,ui) G n x L^(Eo) 



such that, if u is a smooth solution of (3.3) with {u,d^ou) = {uq,ui) on Sq, then 
(3.9) / \u — signjl"^ < Ce, 

and 
(3.10) 



dr Iat' + T^NtqXN' \ . . g 
Here C is a constant that depends on h,F,T but is independent of e. 



elDul'^ + -F(u)] < Ce^. 



7.3 



^Our sign conventions for the unit normal, and hence the mean curvature, are described in Section 
where we also review some basic properties of mean curvature. These sign conventions are such that the curve 
around which the solution in Section 2 concentrates, with the orientation we have implicitly chosen there, in 
fact has "mean curvature" equal to — k rather than k. 



ACCELERATING FRONTS IN SEMILINEAR WAVE EQUATIONS 



11 



The function dr from the theorem is the signed distance to F with respect to the h metric. 
In fact under our hypotheses it is uniformly comparable to the signed Euclidean distance to 



r, so we could replace dr in (3.10) by the Euclidean squared distance with changes only to 



constants (depending on F and the choice of the neighborhood N'). 

Our proof yields additional information that we have not recorded in the statement of the 
theorem, including the following: 

• We show that / elDup + ^F^u) > C > for all small e, so (3.10) implies that the 
energy is strongly concentrated near F. 

• We find certain vector fields X, depending only on the geometry of F, such that 
\\X ■ DuWl^tM N < Ce. 

• Our arguments in fact establish not just estimates of some specific solutions, but also 



more general stability estimates, see Proposition 5.1 for example. 



3.5. discussion. As in [9j and Section [2| the heart of the proof of Theorem 3.2 consists of 
weighted energy estimates in well-chosen coordinates near F. In particular, we will introduce 
coordinates {y^, . . . , y") such that F = {{y^, ■ ■ ■ , y") : = 0}, and in addition y" i— )• c(y") = 
{y^, . . . ,y") is (approximately) a geodesic with respect to the Lorentzian metric for every 
with c(0) € F and c'(0) normal to F. A key point is that the geometry of 



<y 



,n— 1^ 



(y°, . . 

F is exactly such that, when the equation (3.3) is written in the these coordinates, some 



cancellations occur that make very strong energy estimates possible. 



Writing v to denote the solution of (3.3) in the y coordinates, the initial data we consider 
will have the form 



(3.11) 



v{0,y\.. 



<y 



a,ot;(0,y\...,y") = 



near F = {y" = 0}, where q solves 
(3.12) _g" + /o(g) = o, 



/(O) = 0, q{s) — )- ±1 as s — ;> ±oo. 



Existence of a profile q solving (3.12) is standard. Indeed, multiplying by q' , integrating, 



and using (3.4), one finds that the unique solution of 
(3.13) g'_/i(g) = o, g(0) = 



also satisfies ( |3.12 ). We remark for future reference that standard ODE arguments, using 
the fact that /((±1) > (see (3.5)), imply that 



(3.14) 



Iq'(s) - sign(s)| < Ce as s — ^ ±oo. 



As is well-known, the profile q^ is characterized by an optimality property. Indeed, for 
any g : M ^ M, 



by ([3^. Thus \iq{s 
(3.15) 



Co := 



/i(5) 
-7- ±1 as s 

h {s) ds 



- 2^ 



2^ + -Fiq) 



±oo, then 

fliq{s))q'is)ds < 



'-q'^ + lF{q)ds. 



fi (q) ) which occurs exactly when q is a translate 



Moreover, equality holds if and only if eq' 
of the profile qs above. 

Formal arguments suggest that the solution v with initial data (3.11) should satisfy 
v{y^ , . . . ,y'^) ~ leiv"")- As in Remark 2.3 our basic estimate (3.9) in fact implies that 
/[f (y°, . . . ,y'^) — qeiy^)]^ < Ce, but is not sharp enough to distinguish the shape of the 
profile; that is, it does not allow us to say whether sign(y") or q'e(y") is closer to v. But, 
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again as in Remark |2.3[ estimates established in the course of the proof in fact imply that 
for most {y^, . . . , y"'~^), 



dy''' 



for some translation y^ = yQ^y^, • • ■ ^) such that |?/q | < Ce. Indeed, this follows from 
estimates we establish in Proposition 5.1 of ("i, defined in (5.4) below, together with spectral 



estimates like those discussed in Remark 12. 3[ 



3.6. about the nonlinearity. In equation (3.3 ), we have assumed a nonlinearity /e(x, u) := 
fo{u) + £K{x)fi{u), for /o, /i satisfying (3.4), that appears to have a very special form. 

This is not actually the case if fs depends only on u; in this case a nonlinearity fs associated 
to a general double- well potential can be written in this form. Indeed, assume that fs = F^, 
where : M — )• M is a smooth function with nondegenerate local minima at two points, say 
±1, a local maximum at some point in (—1,1), and no other critical points. We claim that 
if these hold, then there exists a number k and a smooth, nonnegative function F satisfying 
(3.5) (and in particular vanishing exactly at x = ±1), such that 

(3.16) /, = F' + eKV2Fsign(_i^i). 

F' and /i 



If we define /o = F and fi = v2Fsign(_]^ then this is exactly (3.4). We sketch a proof: 
For K > 0, let F'^{s) be the unique solution for s > —1 of the ODE (3.16) with initial data 
F'^{—1) = 0. (Despite the non-Lipschitz nonlinearity, uniqueness can be deduced from the 
nondegeneracy of F^, which implies that /^(zbl) > 0.) Then a shooting argument, using 



properties of fs that follow from our assumptions about Fs, shows that there is exactly one 
choice of k such that F'^{s) exists and is positive for s £ (—1, 1), and in addition F'^{1) = 0. 
We then define F*^ outside of [—1,1] by requiring that it solve (3.16) everywhere. Then 
F = F'^ satisfies (3.5) and also solves (3.16) for the value of n found in this argument. 



On the other hand, if fs depends nontrivially on x, then the form of fs has very particular 
useful properties that will be exploited in our analysis. Notably, for fs of this form, the 
optimal profile Qs is formally independent of the value of k. The point is that qs can be 



characterized in terms of either (scaled versions of) (3.12) or (3.13), and hence satisfies 



(3.17) 



e{q 



II , i\ 



1 



fe{q) = -{fo{q) + eKfM) 



for every k S M. (Related issues appear also in (|2.5[), (2.6).) This simplifies our analysis. 



We believe, however, that our arguments could be adapted to study nonlinearities with 
more general dependence on x, at the expense of some technical complications and probably 
weaker estimates. 



4. AN ADAPTED COORDINATE SYSTEM 

We now introduce the coordinate system near T = dl f^ Nt* in which our main estimates 
will take place. 

4.1. a good parametrization of T. Fix a smooth (n — l)-dimensional manifold M diffeo- 
morphic to Fq := F n Sq. Our assumptions imply that M is compact. The example that 
arises most naturally in cosmological settings is M = S*"^^. We will write (y^, . . . ,y"~^), or 
simply y' , to denote local coordinates on M. 
For T > 0, we will write 

(4.1) Mt := {-T,T) X M, 

We will always use standard local coordinates on Mt, by which we mean coordinates of the 
form {y^,y'), where y^ G {—T,T) and y' = {y^, . . . ,y^~^) are local coordinates on M. We 
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will often write y"^ = {y^,y') to denote a point in Mt, where the superscript r stands for 
"tangential" . 

We will parametrize T C Nt* by a smooth map ^ : Mt* — )• Nt* of the form 

(4.2) ^'(y^) = (?/°, V'(y^)) for some ^p : Mt- M". 

(Here and below, we use the fixed coordinate system on to identify it with R^"*"".) Note 
that with this convention and condition (3.2) on the metric {hap)-, assumption (3.7) becomes 

(4.3) ^(0,2/0 = forally'GM. 
We also impose the condition 

(4.4) = = for a = 1, . . . , n - 1, 
where here and in what follows, we use the notation 

(4.5) :=h( — ,^) for a,6 = 0,...,n-l. 

The existence of such a parametrization ^ is rather standard. Indeed, suppose we are given 
local coordinates y' = (y^, . . . ,y"~^) on a subset of M. By definition M is diffeomorphic to 
Fq, so we may fix a diffeomorphism y' — )• ipo{y') S Fq. Then for every y' and sufficiently small 
5 > there exists a unique curve p = p{-; y') : {—5, 6) — )• A^, with components {p^, . . . ,p"), 
such that 

p\t)=t, p{0;y')=My'), Mp'(0,^)=0 for a = 0,...,n-l and |t| <5. 

Indeed, in coordinates this is just a first-order ODE for p(-) to which standard theorems 
apply. It is then easy to see that 5 = 5{y') is bounded away from zero on small enough 
subsets, and on those sets can define ^{y^^y') := p{y^,y')- 

4.2. almost-normal coordinates near the hypersurface. For r positive, we will write 

(4.6) M'' := M X {-r, r). M^ := Mt x (-r, r). 

As above, we will always use standard local coordinates on these spaces, that is, coordinate 
systems that respect the product structure. Thus, in these coordinates, points in M^ have 
the form {y'^,y^) = {y^,y',y^), where y' are local coordinates on M, \y^\ < T and < r. 

The next proposition introduces a map (/> : M^^ — t- A^ that parametrizes a neighbourhood 
of F and such that certain good properties are enjoyed by the pullback metric 

9ai3 ■■= h(^, ^), (s""^) := (5a/3)"\ 9 ■■= det{ga(s) 

where a, /3 G {0, . . . , n}. The proof will be deferred to Section [Tj We remark however that cf) 
essentially defines a Gaussian normal coordinate system near F, modified slightly to arrange 
that condition ( |4.10 ) below holds; this condition implies that changing variables using (p 
maps Cauchy problems for (3.3), with data given for = 0, to Cauchy problems with data 
given on the hypersurface {u^~= 0}. This will be useful. 

Proposition 4.1. For every Tq < T* , there exists (j) : M^^ — )• A^, for some T G {Tq,T*) 
and p > 0, such that such that (j) is a diffeomorphism onto its image, and the following hold. 
First, (j){y'^,0) = ^{y'^), which implies that 

(4.7) TnNT, = HMTo^m, 
and hence that 

(4.8) A^' := 0(M^'') is an open neighbourhood ofTD Ntq ■ 
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In addition, 

(4.9) E Ml'' and = T, then |x°| > Tq for x = (t>{y). 
Second, 

(4.10) X M^/') C So = {x G AT ^ m1+" : x° = 0}. 
Third, the metric satisfies 

(^in block form), where {■fab) was introduced in (4.5). Hence 



(4.12) = ( <^°''(»') °) + 



ly''! (y 



n\2 



In addition, 

(1 
|y"l (y")' 

A''ea;i, i/ie eikonal equation (3.8) /las a unique smooth solution dr on N' , and if we define 
7r"(y'^, . . . , y") = y", i/ien 

(4.14) TT^'o^-^ = dr + 0{dl). 
Finally, 

(4.15) = ^a,.V=y + 0(|y"|) = -^(yO + 0(|y"|). 

Remark 4.2. T/ie implied constants in the above estimates could in principle depend on the 
choice of local coordinates for M . However, our assumptions imply that M is compact, and 
so we may once and for all fix a cover of M by a finite collection of coordinate neighborhoods 
Ui, . . . ,Uk, such that all the above estimates are uniform on (—T,T) x Ui x (—2/9,2/7) for 
every i. (This last point will be evident from the proof of the proposition.) We can then 
require that all subsequent computations in local coordinates are carried out in one of these 
fixed coordinate systems. Having done this, all the above constants are uniform on M^'' . The 
same remark applies below as well. 

Our later energy estimates will contain a symmetric tensor {a"^), defined by 

(4.16) la'^^Up := -g'^'Uo + ^y^'^Ca^/? = -^y^^^o' + ^ E 9''^^^J■ 



It follows from (4.12) that if p is taken to be small enough (which we henceforth assume to 

be the case) then there exists some positive constants 02,03,04, C5 such that 

(4.17) 



^ n— 1 n— 1 

2 E +(l + (y")')e^ < (l + C2(y")')a°^eaC/3 < 2 E ^"'^a^<^ + (1 + C3(y")2)e^, 



a, 6=0 a, 6=0 

n—1 n—l 

(4.18) E ^"'^-^'^ ^ C4 E ^"'^-^b, 

a, 6=0 a, 6=0 

and 

(4.19) ly^^^aCol < |a"^eae/3 

["^P fr.r all f: C ^1+" 



everywhere in MjT, for all ^ G 
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In this section we study the wave equation 



(5.1) 



1 



eU^v + -h{v) + Kfi{v) = 0, v.M^ 



2p 



4.1 



K is a smooth, 



ves (3.3) on N then 



Here g is a metric on M^^ satisfying the conclusions of Proposition 
bounded function on M^'', and /o, /i satisfy (3.4). In particular, if u so. 
V := uo (j) solves (5.1 ) on M^^. 

To state our main estimates, we need some notation. We start by fixing a smooth volume 
form d{vol)o on M, and we extend it to a volume form d{vol) on M^'' by requiring that 

(5.2) d{vol) = dy^ A d{vol)Q A dy" 

in standard local coordinates. We emphasize that d{vol) in general does not coincide with 
the volume form associated to the Lorentzian metric g. 

We will similarly extend d{vol)o to Mt and M'^f, writing d{vol) in every case; the meaning 
should always be clear from the context. 

Thus, in standard local coordinates these are represented by expressions of the form 

d{vol)o = wo(y') dy^ ^ dy"", 
d{vol) = uj{y) dy^ Ady^ A---A dy"^ A dy"" on My'' 

where uj{y'^ ,y' ,y^) = ujoiv') in M^^. Here cjq is a smooth positive function that depends on 
our choice of local coordinates for M. Similar expressions hold for d{vol) on Mt and M^^. 
Next, we define a natural energy density associated to (5.1). For v G H^{M^^), let 

(5.3) eeiv; g) := ^ a'^f'dy.v dy,v + -F{v), 



for F defined in (3.5), and defined in (4.16). We will write simply ee{v) when there is 
no ambiguity, which will be the case throughout this section. Finally, recall that we have 
defined 

Co := y" fi{u) du, 

and that, as noted in (3.15), cq is a lower bound for the energy of a 1-d interface connecting 
the equilibrium states {±1}, and this lower bound is attained by the profile q. 
The following estimates are the heart of the proof of Theorem |3.2[ 



Proposition 5.1. Let v be a smooth solution of (5.1) on M"^^ , and assume that g satisfies 



the conclusions of Proposition \4-l\ Define p{s) = p — c^s, for C5 defined in (4.19), and 



(5.4) 
(5.5) 



Ci(^) : = 
C2(^) : = 



il + C2{ynnee{v)divol) 



CQVoloiM) 



\v - sign{y'^)fd{vol) 



(5.6) C3(^) := / ^ E V V + (y")' Ihdy^M' + -F{v) 



d{vol) 



Then there exists a constant C , independent of v and of e £ (0, 1], such that 

(5.7) Ci{s) < Cmax(Ci(0),C2(0)) fori = 1,2,3 and < s < si := min(r, p/(3c5)) 
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Note that the constants in Proposition 5.1 may depend for example on volo{M), ||k||oo) 
constants in Proposition 4.1 (which may depend on Tq), but they are independent of e G (0, 1]. 

5.1. differential energy inequality. 



Lemma 5.2. Assume the hypotheses of Proposition 5.1. Then 

n-l 



(5.8) -QyO^eiv) < eC( ^ a°''^dyc,v dyi3V + \y''\'^ \dynv\'^) + e dW^ip f + K[evyn - fi{v)] 



VyO 



a,b=0 



g'^'vya ■ Vyo 



where 

(5.9) := (V9^ . . . ,(^"), 99*.- y ^y^ „y„ 

and divj^/2p denotes the divergence on M'^p with respect to the fixed volume form d{vol), so 
that 

(5.10) 

in standard local coordinates on M'^p . 



1 

1=1 



Note that if we compare (5.8) to (2.11), then the term here corresponding to "Term 1" in 



(2.11) has been absorbed into the first term on the right-hand side of (5.8) 
Proof. In standard local coordinates on our equation (|5.ip takes the form 



-A 



-9 



{y^gg^^dy.v) + -Jo{v) + Kf.iv) = 0. 



We rewrite the leading term as a divergence with respect to d{vol) on M^'', leading to 

e 

( 



-dy. (ivg^'^dypv^ - eb'^dy.v + -Jq{v) + Kfi{v) = 



where 
(5.11) 



= g Oy<^ 

-9 



UJ 



Multiply this by f„o and rewrite, recalling (3.4), to find that 



dya{ijjg°'^Vyfi Vyo) + eg^'^Vyfj VyOya + -F{v)yO = eb°'vyaVyO — Kfi{v) VyO 



Note that 



g°'^Vy0VyOy^ = ^dyo{g'^PvyfiVyc.) - ^g^y^VyfiVyc. 



a/3 



and that 9^0 a; = dy^ijj = on . We use these facts and collect all the terms of the form 
dyO [• • • ] to the left-hand side to obtain 



(5.12) 5„o 



-eg^^Vyfi ■ VyO + ^g^'^Vyfi Vy^ + ^F{v) 



^di{ujg'''^dyi3v) + ^gyoVy^Vyp + eb^Vy^VyO - Kfi{v) VyO. 



By definition, the left-hand side is just dyoee{v), and the first term on the right-hand side is 
exactly edivj;/2p (p. It follows from (4.13) and (4.18) that 



ra-l 



gyo Vy-Vy^ <C 2^ a Vya 



a,b=0 
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To estimate the remaining terms on the right-hand side of (5.12) first note that 

n-1 

eb^VyaVyO — K,fl{v)Vy^^ = £ ^ b"'VyaVyO + e{b^ — K)VynVyO + K[£Vy^ — fl{v)]VyO . 

a=0 



Also, since dynoj — 0, we see that ( |4.15p states exactly that — k\ — 0(|?/^|). Thus ( |4.18[ ) 
implies that 

n-1 

1(6" - K)vynVyo\ < C\y'^VyriVyo\ < ^ a°-^VyaVyb + (y")^bj/"P 

a,fe=0 

and that 



n-1 



n-1 



a°'^Vyo.Vyb. 



a=0 



a,b=0 



Thus the Lemma follows by combining these facts with (5.12). 



□ 



5.2. stability of the profile. In this section we collect a couple of lemmas that encode 
some stability properties of initial data for which 0(0) is small, i = 1,2. These concern 
functions of a single variable, which we will denote y", since later we will apply these results 
to the functions y" i— )• v{y'^,y") for y"^ S Mt fixed. 

Lemma 5.3. There exists a constant cq = cq{p) such that if v ^ H^{—p,p) and if 

|y"| |^;-sign(y")|2dy-<c6 

-p 



(5.13) 

then 
(5.14) 



ee,.(t;) dy" > CO - Ce-'=/^ 



for ee,u{v) := -\dyr^v\'^ + -F{v). 



ee,u{v) dy"^ - Co < cy 



Moreover, there exists a constant c-j = c-ji^p) > such that if ( |5.13 ) holds and if 
(5.15) 

then for every p > p, 
rp I 



(5.16) 



'evyn 



fi{v] 



dy^ < C (^J[ee,uiv)dy'' - co) + Ce"^/^ 



as long as v is defined on {—p,p), with C independent of p. 

These are largely proved in Lemma of [9], but since we have modified the statement here 
in some ways, we present some details. 

Proof. Steps 3 and 4 of the proof of Lemma 11 of |9j show that if ce and C7 are chosen to be 
suitably small and (5.13), (5.15) hold, then there exists a function vi and points s_ < s+ in 
(— p, p) such that 

(5.17) |'(;i(s±) - ±1| < Ce-^/^ vi{s) = v{s) \i\s\>p, 

and 



(5.18) 



/p f-p 
ee,u{vi) < / ee,u 
-p ^p 



v). 



In fact, vi is found by minimizing w J Cs^uiw) among the space of functions that agree 
with V outside of the intervals {—p, —jp) U (|p, p). Then (5.18) is clear, and a maximum 



principle argument, together with (5.13^, (5.15 ) and the choices of ce, cy, can be used to show 
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that ±vi{±'^p) > 1 — Ce '^1'^ . If < 1 we can take s+ = and otherwise we can find 

some s+ near \p where = 1. The choice of s_ is similar. 



Let (5(s) := /q* dt^ and note from the definition (3.4) of /i that 
(5.19) e,,,[w) > y/2F{w)\dy^w\ = \fi{w)dynw\ = |9j,nQ(u;(y"))| > OynQiwiy'^)) 

for every w. Thus 



(5.20) / e,,,(w;)(iy" > \Q{w{b)) - Q{wm 

J a 

for every w G i?^ and every a < b. Applying this inequality with w = vi and a = s_, 6 = s+ 



rather easily yields (5.14); see [9] for a little more detail. 



To prove (5.16), assume that v G H^{{—p,p)) for p> p, and note that since F = ^ff 



2 \ V ^^^y" 



eeAv)dy'' - Co + (co - Qlt-Cp)) + (3(^^(-/5))) 



Since ?; = f i at ±p and cq = (5(1) — Q(— 1) = Q{v{s+)) — Q{v{s-)) + 0{e ^/^), we again use 
( |5.20| ) to find that 



CO - QHp)) + QH-p)) < [Q{vi{s+)) - Qiviim + [Q{vi{-p)) - Q{v,isA)] + Ce-^l' 



eeA^i) + / e-sAvi) + Co ) - Co + Ce 

P Js+ 



And again using the choice of s± and (5.20), we have 
Co 



< Q{vi{s+))-Q{vi{sA) + Ce-''' < j^\eAvi) + Ce 



'c/e 



The proof of (5.16) is completed by combining the last three estimates and recalling (5.18). 

□ 

Our next result is exactly Lemma 12 in [2]. Here, in view of future applications, we write 
?; as a function of two variables, G (0)''") ™d y" G {—p,p). 

Lemma 5.4. Let v G H^{{0,t) x {—p,p)) for some r > 0. Then there exists a constant C, 
depending on p hut independent of t and of e £ (0, 1], such that 



|i;(0,2/")-t;(T,2/")|2dy"<C 



(-p,p) 



(0,t)x(-p,p) 



\vl.^'^F{v)dy^dy^. 



5.3. weighted energy estimate. Now we give the 



Proof of Proposition 5.1. We will write Co '■= max (Ci(0), C2(0)) . 
Step 1. Since 

\v{s,y,y^) - sign(y")|2 < 2 {Hs,y,y^) - v{Q,y,y^)\'' + b(0,y,y") - sign(y")|2) , 
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r-p/2 
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C2(5)<2C2(0) + 2 / r \y-\\v{s,y\y^)-v{Q,y',y^fdy^d{vol\ 

JM J-p/2 

< 2Co + C 



p/2 



M \J0 J-p/2 



^-vlo + ^F{v)dy''dyA d{vol)o 



(5.21) <2Co + C [ C3i^)da. 

Jo 



Step 2. For the next few steps of the proof, we regard s as fixed, and we write v{-) rather 
than v{s, •). 

We will say that a point y' G M is good if 
rp/2 

|y"| \v{y',y^)-sign{y^)\'dy^<c,{p/2) 

-p/2 



(5.22) 
and in addition 
(5.23) 



Pis) 

-Pis) 



eeAv){y,yldy^-co<C7{p/2). 



where cg, were fixed in Lemma 5.3 We will say that a point is bad if it is not good. 
We claim that 

(5.24) volo{{y' £ M : y' is bad}) < C(Ci(s) + (2(5)) + Ce"^/". 

To prove this, first note that by Chebyshev's inequality, 

r-p/2 



volo{{y' e M : K2^ fails}) < — / 

= CUs] 



y"| b(y',y")-sign(y-)|2dy"d(7;o0o 



A/ J-p/2 



(5.25) 
Next, note that 



-Ce if y' is good 
/ ,.(f )(s, y\ y'^)dy'^ -cq>\cj if y' is had and (|5.22|) holds 

J-o(s) 



-Co 



always, and in particular if (5.22) fails at y , 



Ci(s) > 



(5.26) 



where we have used Lemma 5.3 for the first case. We integrate to find 
f (I ee,v{v){s,y' ,y'')dy'"' - co\d{vol)Q 

JM \J-p{s) J 

> -Co volo{{y' G M : ([532]) fails}) 

+ C7 volo{{y' G M : y' is bad but ([532]) holds}) - Ce"'^/''. 

Using (5.25), we deduce that 

(z;o/)o({y' G M : y' is bad but holds}) < C (Ci(s) + (2(5)) + Ce'^/', 

and this together with ( 5.25| ) implies (5.24). 

Step 3. Next we estimate C3{^)- We claim that 

(5.27) Csis) < (lis) + CC2is) + Ce-'^' 

for every s. The choice (4.17) of C2 implies that 

n-l 



(1 + C2(y")2)e,(t;) > - J] a'^^Sj^.t; 9^.^; + (l + {yn')eeA^ 

a,b=0 
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By the definitions of (i and ^3 , it follows that 



Ci(^) > ICsis) 



^e,u{v) d{vol) — Co VoIq{M). 



So to complete the proof of (5.27), it suffices to show that 

Co voh{M) - [ ee,u{v) d{vol) < 0(2(3) + Ce-"/', 



(5.28) 



and this follows directly from (5.26) and (5.25). 
Step 4. We now claim that 

(5.29) C'lis) < CiCiis) + C2{s) + C3{s)) + Ce-^/^ 

Recalling the definition p{s) p — C5S, we compute ([(s) = h — C5I2, where 



{s}xMP('') 



l + {y''f)dyoee{v)d{vol) 



{s}xM 



'l + {y^f)ee{v)d{vol)^ 



y"=p{s) 



y"=-p(s) 



It follows directly from the differential energy inequality of Lemma 5.2 that 

(5.30) Ii<CUs) + ha + hb 

where 



ha---- 



{s}xMpM 



lib 



e(l + C2(2/")^) divAfp ip d{vol), 
e(l + C2(?/")2) K (^vyn - ^fi{v)^ d{vol). 



Step 5. To bound Ii^, note that 

(1 + {y^f) divM^P = div,,.p ((1 + (y")2)c^) - 2y>^ 



by (5.10), since a; is independent of y". It is easy to see from the definition (5.9) of ip that 

n-l 



|2/"I \9'"^dyo.v dyOv\ <C a''''dyavdybV + e-^\y''\'^eeA^ 



a,b=0 



and it follows that 



'{s}xMpM 

For the other term, note that d{{s} x M^^^"^) = {s} x M x {—p(s),p{s)} (appropriately 
oriented), and that the induced volume form on dM is just {voI)q. This yields 

y"=p{s) 



{s}xMpM 



divA,2p (1 + (y") V) d{vol) =e (l + (y")^) |(^"| d{vol)o 



{s}xM 



y"-=-p{s) 



Next, the choice (4.19) of C5 was arranged exactly so that ejv'"'! < Cb^e{v), so it follows from 
the above that 



(5.31) 



ha-C^h<CUs). 
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Step 6. We estimate lib as follows. First, 
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lib 



< 



e\\K\ 



{s}xMpM 



Vy0\'^ + [Vyn - ^fl{v) 



d{vol) 



< CUs) + C f (^evy^-^-^\ 



d{vol). 



Note that Lemma 5.3 implies that if y' is good in the sense of (5.22), (5.23), then 



-p{s) 



h{v) 



I 



at y' . Integrating this, we find that 



J{y'eM:y' is good} J-p{s) 
< C 



'evy-n 



Pis) 



dy"^ d{vol)o 
eeAv)dy'' - Co] d{vol)o + Ce-'/' 

' {y'eM:y' is good} W— p(s) / 

< c(Ci{s) + CO volo{{y' e M : y' is bad})) + Ce'^/^ 



Combining these estimates with (5.24), we conclude that Iif, < C(Ci(s)+C2(s)+C3('5)+e~^/^), 
and this together with (5.30) and (5.31) implies (5.29). 

Step 7. Having proved (5.29), (5.21 ), and (5.27), the conclusions of the Proposition follow 
by a Gronwall inequality argument, exactly as in the proof of Proposition 10 in [9j. □ 

6. Proof of Theorem 13.21 
In this section we prove our main theorem. 

6.1. construction of initial data. We will prove that the conclusions of the theorem are 
satisfied by a smooth solution u : Ntq — )• M of (3.3) with initial data 

(u, 9a;0u) 1^0=0 = (^^Oi^^i) constructed below. 

To define uq,ui, we first define (po : M'^p — )• Sq by requiring that 

(/>(0,y',y") = (0,,/.o(2/',2/'^)). 



This definition makes sense in view of (4.10). Next, we set 

(6.1) no = signjg, lii = in Sq \ Image((/)o), 

where Iq := {x £ T,q = M" : (0,x) G /}. In Image((/)o), it is convenient to specify initial data 
in term of the y coordinates introduced in Proposition |4.1[ We would like v = uocp to satisfy 

(6.2) v = vo, dyov = in M'^p, 
when y^ = 0, where 



(6.3) 



vo{y', y") = qeivn ■■= Xpivlqi^) + (1 " Xpiv^) sign(y") 



and xp £ C°°(M) satisfies (2.24). Thus, we complete the definition of uq by 

(6.4) ^0 = ^0° (Po^ Image(^o)- 

We specify ui in Image((/>o) by requiring that 

= dyOV = {dxO'U o (p) dyocj)", 
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when = 0. Thus we define ui = dr^ou\r^o^Q in lmage(0o) by the identity 

n 

(6.5) {ui o (f)) dyOctP = - Y{^^^UQ o (j)) dyO(t>' in {0} x Mp. 

1=1 

The construction of (j) impUes that dyo4P does not vanish in {0} x Mp, so ui is weh-defined, 
and (6.2) holds as desired. 

Observe that the definitions imply that uq and ui are constant, hence smooth, near 
9 (Image ((/)o))- As a result they are smooth everywhere. 



6.2. first estimates of v. As remarked above, v = uocj) solves (5.1) (where we are abusing 
notation and writing k in place of k o (/>) with initial data satisfying (6.2) and coefficients 



of the metric tensor satisfying the conclusions of Proposition 4.1 Thus Proposition 5.1 



is applicable, and in particular v satisfies (5.7). A routine computation, using (3.14) and 
the rightmost inequality in (4.17), shows that for this choice of initial data, the quantity 
max (Ci(0), ^2(0)) appearing in (5.7) satisfies Co < Ce^. In particular it follows that 

n-l 



(6.6) Us) 



a^^dyavdybv + iy 



n\2 



-\dy.v\'' + -F{V) 



for < s < Si = max(T, p/(3c5)). 

6.3. short-time estimates of u in the transition region. Next we define 

e,(n;h) := | i^h'^^d^.uf + h'^d.,rud,,^ + -^F{u). 
Note that our assumptions (|3.1|) on the metric h imply that 



d{vol) < Ce^ 



in M^^. 



-/i°°(a^oM)2 + Y h}^d^^ud^,u ^ \Du\^ := ^(a^.<.n)2 

i,j=l a=0 

in N , where A B means that there exists some constant C such that C^^B < A < CB 
pointwise. One can also check that 

(6.7) ee{u-\i)o(j3^ ee{v;^) 

We next define a smooth cutoff function x" '■ — ^ 1^ such that = 1 outside of 
N' = Image((/)), and in N' , we require that 

o (j){y) = 1 - Xp(y"), where Xp : ^ is defined in ([234]). 

Then is smooth, and x" = near T. In particular, x" ° (kiv) = if \y^\ < g/?, and 

Vx" o 4>{y) = unless \p < < \p- Now fix ti > such that 

(6.8) 

{x E iV : < x° < ti, Vx"(x) / 0} C |(/.(y) : y G M^^ < y° < si, < < . 



Using (6.7) and (6.6), we deduce that 

ee{u;h) dx < C 

<c rc3{s)ds 

Jo 



{x&N:0<xO<ti,Dx^¥'0} 



{yeMl" :0<yO<si ip<|yn|<: 



ee{v;^) d{vol) 



(6.9) 
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Figure 2. 



Here and below, we use Remark 4.2 to obtain uniform bounds on the Jacobian arising in the 
change of variables. 

6.4. energy estimates for u. Computations very similar to those in the proof of Lemma 
15.21 show that 



ee(n;h)— ^ d^i{£h'^d^jud^ou) = -{d^oh°'^)dx<^ud^i3U + eb'^ d^pu d^ou — Kfi{u) d^ou 



where 



2 

1 



-h 



Since fi{u) = 2F{u) by definition, it easily follows that 



(6.10) 



For < r < ti, we deduce that 



C(r) := / xV(n;h) 



l|6.10| 

V C(o) + 



JT. 



<C(0) + C r / x"ee(n,h) +||Dx"l|oo / 



{x:0<2:0<t,|Dx"Ka;)^0} 



<Ce^ + C / C(i) 



In the last line we have used (6.9) and an estimate of C(0) that follows easily from our choice 
of initial data. The integration by parts is justified since u{t, ■) = —1 outside a compact set. 
Hence we can conclude by Gronwall's inequality that 



(6.11) 



(^{t)= [ x"ee(u;h) < Cs^ for all < t < ti. 



6.5. iterate. We now define ■ of the form x" {y^ ly'^y^) = x(|y"l)) where 

XeC-(M), 0<x<l, x = lm(^p,p), supp(x) C (^,2p). 



By arguing as in the proof of (6.9), we deduce from (6.11) find that there exists some 
So G (0, si) such that 



{y(LMl'':0<yO<so,p<\y^\<2p} 



ee{v;g)divol) < Ce^ 
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By combining this with (6.6), we find that 



{j/eA4'':0<i/0<so,DxV0} 



ee{v;g)divol) < C7e2. 



the proof of (6.11), we find that 



Then by using a differential energy inequahty satisfied by v, see (5.12), and arguing as in 

e,{v;g)d{vol) <C£\ 



{so}xMx{y":^<\y^\<p} 



RecalUng from Section 6.2 that Ci('5o) < C'e we deduce that Ci(0; so) < Ce , for 



Ci{s;so):= / il + C2iyn')es{v)divol) 

JmpM 



(6.12) 

We also know from Section [6.21 that 
(6.13) C2(0;so) <Ce2, 



J/°=(s+so) 



covoloiM) 



for C2(s; So) := C2(s - so). 



Now we have shown that v\yO=sO satisfies estimates of the same form (though with larger 
constants) as v\yO^Q. We can thus iterate the above argument to extend estimates first of v, 
then of n, to somewhat longer time intervals. We claim that after piecing together finitely 
many iterations, we can obtain the estimates 



(6.14) / ee{u;h) dx < Ce^ 

Jnto\N' 

(6.15) r ( [ {l + C2{y''f)ee{v;g)d{vol)- CovoloiM)] dy^ < Ce^, 
Jo \Jm^p J 



(6.16) r f ""^'W + (y")' 



a,b=0 



i\dynV\'+-j{v) 



d{vol) < Ce^ 



A proof that finitely many iterations suffice is given in in [9^, proof of Theorem 22] for k = 
and in fiat Minkowski space, but exactly the same proof is valid here. The point is that 
the proof only involves piecing together estimates in the standard and normal coordinate 
systems, and the algorithm for doing so applies equally in this situation, since the (/iq/j) is 
uniformly comparable to the Minkowski metric, see (3.1). 



6.6. conclusion of proof. Now (3.10) follows from (6.14) and (6.16), together with (4.14) 



(6.7) and a change of variables. 



The other conclusion of Theorem 3.2, that is the estimate (2.25) of ||m — sign^ ||^2, follows 



from (6.14), (6.16) and a Poincare inequality. We omit the details, which are just a slightly 



more complicated version of the argument used to deduce (2.25) from (2.17) (via (2.19)) and 



(2.28) in the simple model problem in Section 2 



7. Proof of Proposition 14.11 
In this section we construct a map (f) '■ M"^^ — )■ with the properties summarized in 



Proposition 4.1 
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7.1. construction of 0. To start, for y'^ G Mt we define i'{y'^) G T^(^yT-^T C T^(yT)A^ to be 
the unit normal to F, so that 

(7.1) h(p, i>) = 1, r) = for all r G T^(j,-)r. 

We will fix a sign by requiring that v point "into I" , see below. These conditions uniquely 
determine i>, and our assumptions imply that i— t- v{y'^) is smooth. We next define (j) : 
M^l N, for /9 > 0, by 

(7.2) 0(y^y"):=exp^(,.)(2/V(y-)). 

The condition that P point into / means that (f){y'^, y'^) €z I for all sufficiently small > 0. 
Thus 4> exactly determine Gaussian normal coordinates for near F. 

We will sometimes write {<j)^, . . . , </)") to denote the components of (p in the fixed coordinate 
system on A^. Definition (7.2) states that the components satisfy the system of differential 
equations 



\2ia 



where P", a = 0, 
(7.4) 



+ r 



0, 



, n are the components of v and 



(t>^{y\d) = r{y^), 



dy' 



-(y^O) = i>"(yn, 



1, 



a/3 



are the usual Christoffel symbols. 
Finally, we define 

</.(yO,...,y") :=0(yO-a(y',y"),y\...,y") 
where a : M'^^ — )• M is chosen exactly so that (4.10) holds. Thus, we require that a satisfy 
(7.5) ^0(-a(y',y«),y',y") = in M^''. 

The next lemma implies that the definition of <j) makes sense. 



Lemma 7.1. For p sufficiently small, there exists a : Mp 
addition 

(7.6) |cT(y',y")|<C(y")2. 



satisfying (7.5) and in 



Proof. The definitions (7.2) and (4.2) of cj) and ^ implies that (j)^{y'^ , 0) = y^- Thus it is clear 
that dyO(j)^{y'',0) = 1. Since M is compact, the implicit function theorem thus implies that. 



taking p smaller if necessary, there exists a function a : M'^p 

:ompa 

-(y',0) = 



prove (7.6), it suffices (again using the compactness of M) to show that 



(7.7) 



(j(y',0) = 0, 



such that (|7^ holds. To 

/ that 

for every y' G M. 



The first of these assertions is clear. For the second, we differentiate (7.5) with respect to 
y" and evaluate at a point (y', 0), to find that 



''^ 86^ da da 

-(o,y',o) = %M,^)Z^{y',^) = |^(y',o). 



C^y"^"'"'"^ gyO ' ' QyU ' QyTl 

So in view of (7.3), to complete the proof of (7.7) it suffices to prove that z?0(0,y') = for 
all y' G M. But this follows from noting that at points of the form (0, y'). 







h(z^. 



9y0' 



hoQi^^. 



□ 



26 



BERNARDO GALVAO-SOUSA AND ROBERT L. JERRARD 



Our arguments will imply that (f) is locally invertible in a neighborhood of every point of 
Mt*, and it follows from this that for every T € (To,T*) there exists p > such that 4> is a 
diffeomorphism of Mj!' onto its image. We henceforth assume that this holds. We will also 
feel free to decrease the size of p throughout our argument. 



7.2. estimates of components of the metric tensor. We next prove (4.11), (4.12), and 
(4.13). We will use the notation 

84) d(j) o 



for a, /3 = 0, . . . , n. We first remark that 



g := det{gai3) 



(7.8) 



{lab){v'' 





+ 



0(|y1) 





(in block form), where {^ab) was introduced in (4.5). The estimate 'gabiu^ iV^) = labiv^) + 
0{y'^) for a, 6 < n is immediate when = 0, since (j){y'^,0) = ^{y'^), and then follows by 
the smoothness of (p. The claim above that gem = gn,a = ^na for all (y'^, y") is standard; see 
for example Wald llTj Section 3.3. 

It is convenient to define (t>{y^ ■, ■ ■ ■ , y") := (y" — cr{y' , y^),y^, ■ ■ ■ , y"), so that (p = (p o (p. 
Then the definitions imply that 



(7.9) 



It follows from ^ that g^M{y)) = g^,{y) + 0{{y^f) and that ^ = 0(|y"|) for i > 1 
and from (4.4) and (7.8) that goi = 9i0 = 0{\y"'\) for i > 1. From these one can check that 

9af}{y)=9M + Oi{yn^). 



It follows from this and (7.8) that 

(7.10) {9ap)iy^,yn = 



1 



+ 



Om) 0(|y"|) 
Om) 0{{y^f) 



which is (4.11). Then (4.12) follows by elementary linear algebra considerations; see for 
example Lemma 26 in [9j. 

We must also estimate dog°'^. To do this, we follow ^ and differentiate the identity 
y'^^yp-y = <^7 s-iid rearrange to find that 



Next, after differentiating (7.9) and using (7.8), some calculations show that 
and then (4.13) can be deduced from this and the estimates (4.12) for {g°'^) found above. 



7.3. mean curvature in almost-normal coordinates. We next prove (4.15), which states 
that 

(7.11) ^5-a,.7=9 = -^{yl + o(l?/"l). 



where the left-hand side is evaluated at {y'^,y^) and —K{y'^) is the mean curvature of F at 
(i){y^), since F is assumed to satisfy ( |3.6[ ). Since the left-hand side is smooth, it suffices to 
check this when y" = 0. 
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To do this, we simply verify that the left-hand side of (7.11) equals the mean curvature 
of r with respect to the y coordinates. In these coordinates, F is simply the set {y G M^'' : 
= 0} = Mt, and the unit normal to F by construction is the vector -r^. In general, the 



mean curvature is the trace of the second fundamental form A, where A : TyV x TyT 



IS 



defined by 



A{X,Y) ■.= g{u,VxY). 



Thus 



mean curvature 



^^A{ 



d 



dyo^ ' dyP ' 



«/3 V e — 

dy"- ' dyP ' 



o/3-pn 
9 ^ at 



And from the standard expression for the Christoffel symbols, see (7.4), and using (4.11) 
with = 0, we compute that 



(7.12) 



1 



9^ 



On the other hand, a standard computation shows that 



9al3- 



-1 



-9 



9 OyO 



1 



-9 



„na Ill/a „ 

9 9 ClyO:gfj^i,, 



and by (4.12), this agrees with the right-hand side of (7.12) when y" = 0. 



Remark 7.2. The above proof shows that, in addition to (|7.8|), we have 







1 



+ 



o(|y"l) 





and 



-1 



-9 



-1 



9 



a. 



0(1) 




-9 



9 



<yn+om). 



In particular the first two of these facts, together with (7.8), are somewhat better than the 
corresponding properties of {gap)- 

7.4. solution of the eikonal equation. Finally, it is well-known that if we define dr '■= tt^o 
(f)"^ 1 where ■K^{y^ , . . . ,y") = y", then d^ satisfies (3.8) in Image(<?!»). Indeed, the definitions 



imply that dp = on F, and we have fixed signs such that dr > in Image((/>) n /. Finally, 
by inverting the definition 

and expressing {D(/})~^ in terms of D{(l)^^) o we find that 

o = h''''d,,{r^rd,.{ry. 

Thus the fact that = 1 states exactly that h°'^dx^dr d^/adr = 1. 



: (j) and = Z)(/> everywhere on {y"" = 0}, and it 
follows from the inverse function theorem that = ip^^, D^cf)^^) = D{(f)~^) everywhere on 

0((ip), and since d^ = 



Finally, we deduce from (|7.6|) that 
Hows from the inve: 
F. This implies that 



-l\n 



)", we arrive at (4.14) 
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